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We map certain highly correlated electron systems on lattices with geometrical frustration in 
the motion of added particles or holes to the spatial defect-defect correlations of dimer models in 
different geometries. These models are studied analytically and numerically. We consider different 
coverings for four different lattices: square, honeycomb, triangular, and diamond. In the case of 
hard-core dimer covering, we verify the existed results for the square and triangular lattice and 
obtain new ones for the honeycomb and the diamond lattices while in the case of loop covering we 
obtain new numerical results for all the lattices and use the existing analytical Liouville field theory 
for the case of square lattice. The results show power-law correlations for the square and honeycomb 
lattice, while exponential decay with distance is found for the triangular and exponential decay with 
the inverse distance on the diamond lattice. We relate this fact with the lack of bipartiteness of 
the triangular lattice and in the latter case with the three-dimensionality of the diamond. The 
connection of our findings to the problem of fractionalized charge in such lattices is pointed out. 
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I. INTRODUCTION 

There has been enormous interest in the properties of 
quantum frustrated magnets, with recent studies show- 
ing fractionaHzation of quantum numbers (see Refs. 00|, 
and citations therein). More recently, the attention 
has shifted to charge degrees of freedom on lattices 
with frustrated geometries where the number of classi- 
cal ground state configurations increases exponentially 
with the number of sites. It has been noticed that hop- 
ping of fermions from site to site in frustrated lattices 
can lead, under certain conditions, to fractionaHzation of 
the fermionic charge in a rather natural way. -^ 

The motion of spinless fermions (or hardcore bosons) 
on a lattice can be mapped to an Ising spin problem or 
to a dimer configuration for any given lattice geometry 
and filling fraction as we will see in the next paragraph. 
The addition of a fermion (or boson) leads to a fiip of an 
Ising spin or, equivalently, to a new dimer. The ground 
states of the considered systems fulfill a local constraint 
of having a certain number of particles on each unit cell 
or, equivalentely, a certain number of dimers at each site. 

For a specific implementation of local constraints in 
the Hamiltonian language for fermions, let us consider 
the model Hamiltonian 



(ij) 



<«J> 



Here (i j) denotes the sum over nearest neighbors i and 
j. The operators Cj(q) annihilate (create) a particle on 
site i, and rii — c\c^. For a moment, let us set the hop- 
ping integral t to zero. We are interested in lattices for 
which the ground state of the Hamiltonian Q with re- 
pulsive nearest-neighbor interaction term V only has at 



certain fillings a macroscopic degeneracy, which increases 
exponentially with the system size. In other words, the 
system has a finite T = entropy (for possible technical 
applications, see Ref. jSJ). Examples of such lattice are 
the pyrochlore lattice and the criss-crossed checkerboard 
lattice (two dimensional pyrochlore) , shown as Fig. ^ (a) 
and (b) . At half filling, all classical ground states fulfill 
the so-called tetrahedron rule of having exactly two par- 
ticles on each tetrahedron (criss-crossed square )2i. The 
dimer coverings of the dual lattice is constructed as usual 
by connecting all centers of the tetrahedra and drawing 
a dimer whenever the traversed site of the original lat- 
tice is occupied. The tetrahedron rule translates into the 
constraint of having exactly two dimers at each site of 
the dual lattice. Similarly, for the kagome lattice with 
nearest-neighbor repulsion the dual lattice is the honey- 
comb, see Fig. Q](c). If the repulsive interaction on the 
kagome lattice is restricted to hexagons only, the dual 
lattice is a triangular lattice, see Fig.Q](d). 

In the long-run, one wants to understand the behavior 
of mobile, fractionally charged, excitations resulting from 
weak particle or hole doping of such systems. 

As a first step, we present the study of classical defect- 
defect correlations. A defect here means a cell where the 
rule of fixed number of particles (or dimers) is violated 
(Fig. 2). The study includes bipartite and non-bipartite 
two-dimensional lattices (checkerboard, honeycomb and 
triangular lattices) as well as a three-dimensional dia- 
mond lattice, with two different fillings corresponding 
to (i) the loop dimer and the (ii) hard-core dimer cover- 
ing. Wherever possible, we obtain analytical results and 
compare them to numerical simulations. These classical 
correlations provide informations about correlations near 
the Rokshar-Kivelson point (RK point) of the quantum 
Hamiltonian because at the RK point the quantum me- 




FIG. 1: (color online) Mapping of lattices which show frus- 
tration (black) to dual lattices (red), (a) Checkerboard lat- 
tice with nearest-neighbor repulsion -^ square lattice, (c) py- 
rochlore lattice with nearest-neighbor repulsion -^ diamond 
lattice (c) kagome with nearest-neighbor repulsion -^ honey- 
comb lattice and (d) kagome with repulsion interaction on 
hexagons only -^ triangular lattice. The constraint of having 
a certain number of particles on each unit cell translates into 
the constraint of having a fixed number of dimers touching 
each site of the dual lattice, leading to a hard-core dimer or 
loop dimer covering. 




FIG. 2: Hard-core dimer (a) and loop dimer (b) covering on 
the square lattice. Figures (c) and (d) show two defects in 
the different coverings at postions marked by arrows. 



chanical ground state is given by an equal weighted super- 
position of all configurations^. Furthermore, in our case, 
one can expect classical correlations to hold more gen- 
erally since the low-energy excitations can be described 
equivalently by hardcore bosons or fermionsSi. 

Historically, studies of monomer and dimer models 
started very earlj*^. At zero doping, the PfafRan method 
introduced by Kasteleyn allows to obtain classical corre- 



lations analyticalljiiS. For the square lattice with hard- 
core dimer covering, the correlation of a pair of monomers 
has been obtainedii^. 

The main findings of the present work are: For the 
square and the honeycomb lattice, correlations as func- 
tion of distance show a power-law behavior, with an ex- 
ponent 1/2 for the hard-core dimer covering on the square 
lattice and 1/3 for a loop dimer covering. The honeycomb 
lattice has an exponent of 1/2 for both coverings. An ex- 
ponential decay with distance is found for the triangular 
lattice and the diamond lattice shows exponential decay 
with the inverse distance. 

In the next section, we present the results and the 
methods in detail. We start with the square lattice where, 
for completeness, we summarize the Liouville theory of 
the fully packed loop model with two colorsiSi, which pro- 
vides the exponent for the loop covering analytically and 
we subsequently show the numerical results. We then 
present the relevant results for the honeycomb, the trian- 
gular and the diamond lattices and finally, we conclude. 



II. RESULTS FOR DIFFERENT LATTICES 

A. Square lattice 

Square lattice. As described above, the half-filled 
checkerboard lattice with nearest-neighbor repulsion is 
mapped to a loop dimer model on the square lattice, 
meaning that each site is connected to exactly two 
dimers. This can then be interpreted as covering in the 
fully packed loop model with two colors, assigning dif- 
ferent colors to occupied and unoccupied bonds. Using 
the machinery developed for this model, ^^ we can read 
off the appropriate exponent of the defect-defect corre- 
lation function: In order to obtain the Liouville theory 
of the model, the first step is to map the oriented loops 
to an interface model, where they can be interpreted as 
contour lines. The microscopic heights are defined at 
the center of the lattice plaquettes and each bond is in 
one of four possible states {A, B, C, D). We write the 

where Nh and 
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partition function as Z — Yl 
Ng are the numbers of the loops of the two colors and 
< nfc, Ug <2 are the fugacities accordingly. Assigning 
weights exp(±i7ref,) or exp(±i7reg) to the different orien- 
tations, with nil = 2cos(7re6) and rig = 2cos(7reg), allows 
a local definition of the statistical weight of oriented loop, 
and in turn the weights leads to a local field theoryji^ 

Our problem at hand (which can be equivalently 
mapped to a six vertex model), corresponds to fugac- 
ities n&=ng=l. Therefore eb=eg=l/3. Adopting the 
same convention as in Ref. [i3| [A = (— 1,+1,+1), 
B = (+1, +1, -1), C = (-1, -1, -1),^^ =J+1, -1, +1)], 
the topological charge (±m) \s rn — C — A— [Q, —2, —2) 
and 2xi 1 is the dimension of an operator with total 
charge (eoi'Tii^i), where eg = -7r/2(eg -I- ef,,0,eg — Cfc) = 
7r(l/3,0,0). The effective field theory for the coarse- 



grained heights is given by the action 
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Here, Se is the elastic term (entropy of the oriented 
loops) with the tensor K, of elastic constants Kn = ^33 
and K12 = K23 = and h{x) is the height field (coarse 
grained field of the integer- valued height h). Sb is the 
boundary term with the scalar curvature TZ, which van- 
ishes everywhere except at the boundary. This term in- 
serts the vertex operator at far ends and supplies wind- 
ing loops. Finally, Sl is the Liouville term. The Li- 
ouville potential is the coarse grained version of the lo- 
cal weight assigned to an oriented loop configuration 
J|-A(x): u){x) = — ln(A(x)) is defined in height space, 
which is the bcc lattice Ru,. The Fourier-transformed 
values oi U!{x) = X^eGfl* ^^ exp(ie • h{x)) in the recipro- 
cal lattice i?* (a fee lattice) enter Sl- 

The scaHng dimension of a general operator with both 
electric and magnetic charge is given by: 
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where in our case Kn — ^{2 — et 
f (eb - eg) = and i^s" - - ''~'^^^'" 



the critical exponent of interest is simply: 



n ■ [/C • m] (2) 
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The result (PJ will be next verified numerically by Monte 
Carlo simulations. 

The dimer models on dual lattices are numerically 
preferable over the original lattices, because the con- 
straints are included in a more natural way. Even though 
dimer models have been intensively investigated in the 
study of spin model^Si not much is known about defect- 
defect correlations. The defect-defect correlation func- 
tion, e.g. C(0,r) r^ {NoNr) which gives the probability 
to find two defects at a certain distance r. The defect- 
defect correlation is proportional to the restricted parti- 
tion function Zii(0,r) ~ |j.|-2a;ii ^j ^j^g ^^^ color fully 
packed loop model. It is counting the number of con- 
figurations with defects at and r connected by one 
string of each color. Thus C(0,r) is expected to become 
isotropic at large distances and to decay algebraically 
Z(0,r)~|r|i/3. 



In the following, we measure numerically the correla- 
tions along a coordinate axis [r = {x, 0)] and refer to it 
as C{x). The classical two-point correlation functions at 
zero temperature are delivered as averages over all degen- 
erate ground states. The number of degenerate ground 
states grows exponentially with the system size. At the 
start, an allowed configuration with fixed fiUing, i.e. a 
ground state of the undoped system with no violation of 
any local constraint is generated. Then we add a dimer 
onto an unoccupied random link. This leads to two de- 
fects on adjacent sites of the dual lattice, which subse- 
quently propagate via local dimer moves through the sys- 
tem without creating any new defects. At each step the 
defect-defect distance x is counted in a histogram Z{x), 
which after normalization yields the correlation function 
C{x) = Z{x)/Zq. The normalization is somewhat arbi- 
trary, we chose it in such way that C(l) = 1. The algo- 
rithm is terminated when the standard deviation of the 
measured quantity falls below a certain threshold. Re- 
sults of simulations with different initial configurations 
are compared for verification. For the 2D lattices about 
10^ — 10* sampHngs were necessary. For the 3D diamond 
lattice it took about 10^" steps until convergence. We 
also applied an alternative Monte Carlo algorithm with 
loop updateaii which is known to be ergodic and unbi- 
ased but shows considerable slower convergence. 

The algorithms have been applied to different lattice 
structures and fiUing factors (see Fig. QJ. To test the 
implementation of the algorithm, we first reproduced the 
known monomer two-point correlations on the hard-core 
dimer covering on the square and triangular lattice (see 
Fig. and Fig. El as well as the dipolar correlations in 
the undoped system on a square lattice i ^"''^'? 

The square lattice is a bipartite lattice and the two 
defects are on different sublattices. We extracted the 
exponent from the numerical data by linear interpola- 
tion of log-log plots and verified the results by finite size 
scaling C{x/ L) = U' [x / L)'' c{x / L) with exponent 7 and 
system size L. Fig. compares the numerical data and 
the analytical results. For a rescaled distance variable 
X = Lsm{Trx/L)/'!r is used the numerical fit in order to 
account for the periodic boundary conditions. 

In the case of a hard-core dimer covering (quarter-filled 
checkerboard lattice) the exponent agrees with the results 
from Ref. [la|, i.e., C{x) ~ 1/x^'^. The correlations 
in the case of a loop dimer covering (half-filled checker- 
board) are very well fitted with the power law l/x^^'^, 
obtained analytically (see above). As expected, the cor- 
relations do not show any angle dependence at large dis- 
tances. In both cases, the decay of the correlation func- 
tion is algebraic which is expected for two-dimensional 
bipartite lattices: The two defects have long-range corre- 
lation and feel each others presence at all distances. The 
bipartiteness is also seen by the strictly zero correlations 
C(x) for distances which connect sites that reside in the 
same sublattice. 
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FIG. 3: Classical defect-defect correlation functions on a 
square lattice with hard-core dimer (upper panel) and loop 
dimer coverings (lower panel) along a coodinate axis. Results 
are shown for lattices with L = 64 and L = 128. The numeri- 
cal data is fitted by the exact asymptotic results. The periodic 
boundary conditions are taken into account by plotting the 
fit against a rescaled distance variable. 
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FIG. 5: Classical defect-defect correlation functions on a tri- 
angular lattice with hard-core dimer (upper panel) and loop 
dimer coverings (lower panel) along a coodinate axis. Results 
are shown for a lattice with L = 16. 



C. Triangular lattice 
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FIG. 4: Classical defect-defect correlation functions on a hon- 
eycomb lattice with hard-core dimer and loop dimer coverings 
along an axis of the bravais lattice. Numerical results are 
shown for lattices with L = 64 and L = 128 together with 
to a power law ^ x ' . The periodic boundary conditions 
are taken into account by plotting the fit against a rescaled 
distance variable. 



B. Honeycomb lattice 

The kagome lattice with nearest neighbor repulsion 
at one third and two third filling are mapped to the 
hard-core dimer and the loop dimer model on the hon- 
eycomb lattice, respectively. The numerically obtained 
correlations along the axes of the bravais lattice are both 
C{x) ~ 1/x^/^ (see Fig. E||. Note that the two models, 
i.e., filling factors, are equivalent in the absence of de- 
fects. They can be identified by exchanging links which 
are occupied by a dimer and those which are not occu- 
pied. The faster algorithm described above did not show 
the equivalence of the two models. Therefore, we used 
the loop algorithm for this lattice. 



The kagome lattice with repulsion on the hexagons 
only is mapped at one sixth filling to a hard-core dimer 
model and at one third filling to a loop dimer model 
on the triangular lattice. The correlation function de- 
cays exponentially in both cases (see Fig. EJ with decay 
length of the order of one lattice spacing. The values ob- 
tained for the hard-core dimer model are in agreement 
with those obtained in Ref. [iJl- For distances fur- 
ther than a few lattice spacing C{x) is constant within 
the noise ratio and tends to a finite value in the limit 
X — > ooi£. This impHes that the free energy difference 
AF(cx)) T(lnC(oo) - InC(l)) of two infinitely sep- 
arated defects in the classical hard-core dimer and loop 
dimer model is finite and the two defects are deconfined. 



D. Diamond lattice 

The mapping of particles with local constraints due 
to large repulsive interaction to hard-core dimer or loop 
loop models can be applied in three dimensional systems 
as well. For example, the pyrochlore lattice at one fourth 
and one half filling corresponds to the hard-core dimer 
and the loop dimer model on the diamond lattice, re- 
spectively. 

Monte Carlo simulations for two different lattice sizes 
are shown in Fig. |£| Much more samplings are needed to 
achieve a good signal to noise ratio, because the phase 
space is considerable larger than in the case of the 2D 
lattice. From a logarithmic plot, we found that the corre- 
lation functions decay approximately exponentially with 
respect to the inverse distance along the axes of the bra- 
vais lattice as C{x) ~ exp(l/4a:) for the hard-core dimer 
covering and C{x) ~ exp(l/62:) for the loop covering. 
The reason that this bipartite lattice shows these short 
range correlations is because the monomers follow a 3D 
Coulomb law with potential (y{r) ^^ 1/r) as opposed to 



+ /x Dimer, L=32/64 
V/A Loop, L=32/64 




X >"0<A axXaxXx: : 



1 

0.95 

O 0.9 

0.85 

0.8 

5 10 15 20 25 30 

Distance along axis x 

FIG. 6: Classical defect-defect correlation functions on a di- 
amond lattice with hard-core dimer (upper panel) and loop 
coverings (lower panel) along a coodinate axis. Results are 
shown for lattices with L = 32 and L = 64. The numer- 
ical data are fitted to exponential functions. The periodic 
boundary conditions are taken into account by plotting the 
fit against a rescaled distance variable. 



the divergence at large distances V{r) oc logr found in 
2D2fi. The restricted partition function then decays ex- 
ponentially with the distance in the case of the diamond 
lattice. A recent work on spin 1/2 Heisenberg antifer- 
romagnet on the pyrochlore latticeSl reveals a fraction- 
alized spin Hquid with U(l) gauge structure where this 
l/r potential acts between pairs of spinous and pairs of 
monopoles. 



III. DISCUSSION 

In summary, we calculated the classical defect-defect 
correlations for different fillings and different lattices with 
geometry that causes frustration in the motion of par- 
ticles. Recent studies on bipartite and non-bipartite 
three-dimensional lattices^ lead to the conjecture that 
extended critical phases are reaHzed only in bipartite 
lattices. This is supported by our results: For two- 
dimensional bipartite lattice/dual lattice pairs (checker- 
board/square lattice and kagome/honeycomb lattice) the 
correlation functions decay with distance following a 
power law and tend to zero for x — > oo. Other- 
wise ("modified kagome with repulsive interaction on 



hexagons only /triangular lattice and three dimensional 
pyrochlore/diamond lattice") an exponential decay with 
distance or inverse distance, respectively, and a finite 
value for a; ^ c» is found. 

Where analytical values for power-law exponents were 
known in the literature, these could be confirmed. For 
the checkerboard lattice, i.e., two-dimensional pyrochlore 
lattice, we could relate the defect-defect correlations to 
the solved two-color fully-packed loop modeU^. This pre- 
dicts C{x) ~ x~^/^ which is in perfect agreement with 
the data from our Monte Carlo simulation. 



lattice 


hard-core 


loop 


square 


a;-i/2 


a;-i/3 


triangle 


exp{—x) 


exp(—x) 


honeycombe 


x-'/^ 


x-^/^ 


diamond 


exp{l/4:x) 


exp(l/Qx) 



TABLE I: The summary of the defect-defect correlations we 
obtained for the four different lattices and the two distinct 
dimer coverings. In the case of the triangular lattice there 
is an exponential decay with decay distance of one lattice 
constant, tending to a finite value at x — > oo. 



The different behavior of the defect-defect correlations 
leads to two different scenarios with respect to the separa- 
tion of defects at small finite temperatures. A separation 
of two defects to an infinite distance leads to an increase 

of the free energy, Ai^(cx)) T{hiC[oo) - InC(l)). 

Consequentely, AF(oo) is infinite for the defect on the 
square and honeycomb lattice (confinement) and remains 
finite for the triangular and diamond lattice (deconfine- 
ment) . 

The results are useful for the quantum version of the 
problem at the RK point of the Hamiltonian in general^. 
In particular, our results can be extended to the motion of 
strongly correlated spinless fermions on the checkerboard 
lattice^.. 
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